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The proximity-effect theory developed by Takahashi and Tachiki§ for infinite multilayers 
is applied to multilayer systems with a finite number of layers in the growth direction. The 
purpose is to investigate why previous applications to infinite multilayers fail to describe 
the measured data satisfactorilya. Surface superconductivity may appear, depending on the 
thickness of the covering normal metallic N layers on both the top and the bottom. The 
parameters used are characteristic for V/Ag and Nb/Pd systems. The nucleation process is 
studied as a function of the system parameters. 
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An alternating sequence of superconductor and normal metal layers (S/N) generates a system whose properties 
have raised both theoretical and experimental interests. The dirty-limit theory of Takahashi and Tachikld is meant to 
describe such systems. However, it has some insufficiencies in explaining experimental results such as the dimensional 
crossover in H C 2 ,\\(T) curves,au displaying the upper critical field applied parallel to the layers as a function of the 
temperature. Since up to now exact calculations using this theory were done for infinite periodic systems, we extend 
the calculations to finite multilayers. A marked difference is found in the proximity effect on the pair amplitude for 
the upper critical magnetic fields H C 2,j_ and n. In trying to cover as many aspects as possible which underlie the 
physics in such systems, we are hoping to prepare a theoretical explanation of the experimentally observed dimensional 
Q ■ crossover in the parallel upper critical field -ff C 2,|| versus temperature, which improves upon the at present available 
O \ explanation, being a rather artificial ad hoc explanations We apply the theory to real systems such as finite multilayers 
of V/Ag. 

Section || gives a brief survey of the Takahashi- Tachiki theory. In Section III we calculate the spatial distribution 



^ i of the pair amplitude, showing features of the proximity effect, which are typical for layered systems. Two distinct 
' situations are treated separately, namely one in which the magnetic field applied to the system is perpendicular to the 
layers and the other one when the field is parallel to the layers. Section [V is devoted to nucleation properties, which 
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form an intriguing particularity of superconductivity when a parallel magnetic field is applied. The manifestation of 
surface superconductivity and related problems, such as the influence and role played by the boundary conditions, are 
treated. In Section ^ we will present a detailed description of the dependence of the proximity effect on the parameters 



of the system. The results obtained restore confidence in the applied microscopic theory for dirty multilayers, 
-t— > ' 

II. THE TAKAHASHI-TACHIKI THEORY 

*rt ■ . 

In inhomogeneous systems such as S/N multilayers, the superconducting properties are influenced by the proximity 
of different materials. One can think of non-superconducting materials, or of superconducting materials with a lower 
T c . A finite Cooper pair amplitude is induced in the less- or non-superconducting material. This proximity effect 
L| ' leads to a reduction of the critical temperature T c compared to that of a bulk superconductor. 
.,— i \ We will apply the description developed by .Takahashi and Tachiki, which has been devised to be valid for inho- 
mogeneous superconductors in the dirty limitJil We now give a brief summary of the main equations, by that also 
5h ' defining the system parameters to be used. _ 

The theory starts from the Gor'kov equatioro for the pair potential A(r) 

A(r) = V»fcT^ f d 3 r'Q w (r, r')A(r'), (1) 

in which the material parameter V(r) is the space dependent BCS coupling constant and the integration kernel Q u 
obeys a Green's function like equation 

[2\w\ + L(V)]Q w (r, r') = 27rV(r)<5(r - r'). (2) 
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The differential operator is defined by 



L(V) = -fcD(r)(V - ^A(r)) 2 . 



(3) 



The material parameters JV(r) and -D( r ) are the electronic density of states at the Fermi energy and the diffusion 
coefficient, respectively. V(r), N(r) and D(r) are constant in each single layer. At the interfaces de Gennes boundary 
conditions are imposed,i which require the continuity of and D(r)(\7 — ^A(r))F (r), where the pair amplitude 
F(r) is related to the gap function A(r) through 



A(r) = V(r)F(r). 



(4) 



Takahashi and Tachiki provide a way of solving Eqs. ([j]) and (0) by developing the kernel Q w (r, r') and the pair 
function F(r) in terms of a complete set of eigenf unctions of the operator L(V). These eigenf unctions are labeled by 
the parameter A, and the eigenvalues are ex- They are solution of the eigenvalue problem 



L(V)* A = ex^x- 

The requirement of the existence of a solution for Eq. ([!]) leads to the equation 



det|<5 AV -27rfcT^ 



1 



2U +e x 



V\\i | = 0. 



(5) 



(6) 



For finite multilayers covered by infinite N layers, the boundary conditions at infinity, imposed in determining the 
eigenf unctions ^x(x,y,z) of the differential operator L(V), are such that ^\(z — > ±oo) = 0, which is in line with 
F(z — » ±oo) = 0. As usual for these type of layered systems, the growth direction coincides with the z direction. For 
finite multilayers in vacuum, we apply de Gennes boundary conditions, since they insure that there is no current flow 
through the interface between the multilayer and the vacuum. These boundary conditions read ^£id£igig) = at the 
interface with the vacuum. In the absence of a magnetic field, the solution of Eq. (0) giving the largest value for 
the critical temperature is the physical one. In the presence of a field, solving this equation allows us to derive the 
H C 2(T) curves. The temperature at which H C 2 — > is T c . 



III. THE PERPENDICULAR AND PARALLEL UPPER CRITICAL MAGNETIC FIELDS 



We consider a multilayer system such as the one depicted in FIG. [I]. A certain number of finite-thickness S/N layers 
are covered by infinitely thick external iV-layers in the growth direction z. We will vary the number of internal layers, 
and show results for systems with 1, 5, 9, and 69 internal layers respectively. A magnetic field of magnitude 1 kGauss 
is applied in the growth direction, so perpendicular to the layers. The computation and analysis of the pair function 
will lead us to some interesting conclusions related to the microscopic aspects of the proximity effect. 




FIG. 1. The movement of the electrons in the system, when a perpendicular magnetic field is applied. 



For a perpendicular orientation of the field we can choose A(r) = (0, Hx,0). By that, applying the method of 
separation of variables, we can write the solution ^x of Eq. (||) as 

*A=^(x)e lfc ^$ A (z) (7) 

The solutions in the x direction are parabolic cylinder functions, which for the present case are equivalent to one- 
dimensional harmonic oscillator functions. The system is infinite and uniform in the x and y directions, so that the 
ground state has a x — — \ and arbitrary k y . With ^x given by Eq. (^), Eq. (^) reduces to 
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with 



in which the magnetic coherence length £ 



dz 2 



<f x (z) = kl(z)^x(z), 



TiD(z) 
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(8) 
(9) 




FIG. 2. The pair function in perpendicular magnetic field as a function of the coordinate z in the growth direction, 
thickness of both N and S layers are ds = djv = 100A. The magnetic field applied is H±=l kGauss. 
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In FIG. ^ the pair function is shown for four -systems with a different number of internal layers. The system 
parameters used are specific for V/Ag multilayeraJa. Two features of the proximity effect come out clearly. One is 
that the proximity of the S and N layers leads to the oscillating behavior, in which the minima correspond to the N 
layers, while the maxima occur within the S layers. The other one is an overall proximity effect due to the infinite 
exterior iV-layers. The 69-layers result can be considered as to stand for the limit of a large number of layers. The 
influence of the outer layers is not noticed any more in the middle of the system, the behavior here becoming similar 
to that of an infinite periodic system, described by Koperdraad et al.B 
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FIG. 3. The movement of the electrons in the system, when a parallel magnetic field is applied. 

Turning to the parallel field case, we first show in FIG. || in an intuitive way how the electrons move. This leads to 
a movement of the Cooper pairs in which they would extend naturally in the normal layers. But as they are forbidden 
there, the concentration of Cooper pairs in the superconducting layers will be pushed up. As in the perpendicular 
field case, the outer N layers extend to infinity. 

To compute the pair functions, we go back to Eq. ([5]). A parallel magnetic field lies in the a;y-plane, and we can 
choose the vector potential A(r) as (Hz, 0,0). Separation of the variables leads to the form 



*a = e 



Jk y y 



W\(z) 



for \, and the equation for w\ can be written as 

, d 2 (z 
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-W\(z). 



(10) 



(11) 



The point zq = (, k x is called the nucleation point. For the results to be presented in this section the nucleation point 
is chosen in the middle of the system, so zq = 0. In the following section it will become clear, that this symmetry- 
inspired choice is not the only possibility, particularly if the covering N layers are no longer infinitely thick and the 
system becomes really finite in the growth direction. 
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FIG. 4. The pair function for systems with a different number of layers in a parallel magnetic field, for the layer thickness 
ds = djv = 100A. The applied magnetic field is H\\ = 7 kGauss. 

In FIG. | we show the pair functions of the systems containing 1, 5, 9, and 13 internal layers. The markedly 
different behaviour compared to what was found for the perpendicular field case is certainly due to the anisotropy of 
the system. For each system considered, the pair function drops to zero, except for a few layers, centered around the 
nucleation point zq. Apparently now the proximity effect, here being equivalent to the spreading of the pair amplitude, 
is much weaker than for H c2 ±. The result shown is obtained for a magnetic field of 7 kGauss, which is somewhat 
larger than the field used in the perpendicular field case. In the latter case the features of the pair distribution are 
insensitive to the magnitude of the magnetic field. But it appears that in the parallel field case, the higher the field is, 
the more the pair function is concentrated at the center. So one expects not much difference in the H c2 \\(T) curves 
for lower temperatures. This is ilustrated nicely by FIG. ^, in which these curves are shown for the monolayer and 
5-layers systems. The H c2y ±(T) curves are shown as well. The H c2 ^(T) curves indeed coincide for higher parallel 
magnetic fields, irrespective to the number of internal layers considered, but they are clearly different for a small 
parallel magnetic field. This is related to the fact, that for small magnetic fields, the pair function is spread over the 
multilayer, in a similar way as it is for the perpendicular magnetic field. 
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FIG. 5. #c2,||(T) and H c2 ,±{T) curves. 

The fact that the parallel upper critical field is larger than the perpendicular field is a consequence of the difference 
in the distribution of the pair amplitude for the two field directions. In the parallel field case, the pairs are concentrated 
in the S layers. That is why for this direction a larger magnetic field is needed to destroy the superconductivity. A 
similar argument explains the fact, that H c2 ± for five layers is larger than for one layer. Since for ff C 2,_i_ the pair 
amplitude is spread over the layers, a 5-layers system is a better superconductor than a one-layer system for all 
temperatures. Consequently, a one-layer system is expected to have a lower T c than a 5-layers system. This is shown 
clearly by the H c2 j_(T) curves, and that is why the parallel curves converge to the perpendicular ones in the zero-field 
limit. 

Now we turn to really finite systems, for which the covering N layers are finite as well. 
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IV. SEARCH FOR THE NUCLEATION POINT 



Measurements of the parallel upper critical magnetic field H C 2.\\ for multilayer systems such as V/AgtjEl show 
interesting features, such as the dimensional crossover in H C 2 : \\(T) curves. A generic H C 2 t \\(T) curve is divided in three 
regions. Near T c , there is an average 3D behavior, corresponding to a linear dependence of the critical magnetic field 
on the temperature. At a lower temperature, at which the magnetic coherence length becomes of the order of the 
layer thickness, this behavior is replaced with a square- root-like dependence of -ff C 2,|| on T, corresponding to a 2D 
domain. A second dimensional crossover is only present in S/S' systems, which comes-.fapm a difference in diffusion 
constants for the two superconductors, and to which we do not pay further attentionJiiM 

Up to now, the occurrence of the dimensional crossover has not been described satisfactorily. The calculated 
crossover temperature appears to lie much higher than the measured one.0 But the results obtained so far only apply 
for infinite periodic multilayers. Such model systems have translational symmetry in the growth direction, so that 
surface effects are excluded. It was AartSj-who suggested to investigate finite size effects, which always may show up 
in the finite samples used in experiments.113 

An extension of the applications of the theory to finite systems is not a trivial one. The reason is, that in calculating 
the parallel upper critical magnetic field one has to consider the position of the nucleation point zq. In infinite periodic 
systems, nucleation always occurs in the middle of a layer, due to the periodical symmetry of the system. Each center 
of a layer is a mirror-symmetry center. One only has to figure out which layer of the two types of layers will 
correspond to the optimal nucleation. In practice this has led to a distinction of two possible solutions, arfirst solution 
corresponding to nucleation in the S layers, and a second solution describing nucleation in the N layersE 

For finite multilayers the symmetry is largely reduced. Due to mirror-symmetry with respect to the center of the 
system, one would be inclined to implement the choice of the nucleation point in the center only, as we did in the 
previous section. But finite systems are known to exhibit surface superconductivity as well. This means that one has 
to search for the precise location of the nucleation point. If it would turn out to lie close to the surface, it would 
lie at a point devoid of any symmetry. While surface superconductivity is well investigated in finite homogeneous 
superconductors, as far as the authors know, such a search has never been done for inhomogeneous layered systems. 
This is the subject of the present section. 

In looking for the correct solution of Eq. (||) it is important to make use of the fact, that the eigenfunctions ( |l0| ) 
corresponding to different nucleation points are orthogonal, since they have different wave-numbers k x . Due to this, 
they provide zero-matrix elements in the matrix V\\>. The determinant will be split in an infinite number of 
blocks, each block corresponding to a certain nucleation point zq = £ k x . 

det \5 XX , - 2nkTj2 9| \ V xx ,\ = Rdet \S x , Zo ., x ,, Zo - 2nkTj2 ^JT Vx^jV^I = 0- (12) 

— i \ui \ + e x — i \U! \ + e x .z a 

UJ 1 1 20 UJ ' 

In the determinant equation (|l^), we are interested to solve for that block 

det \S XtZQ , x ^ Z0 - 27rfcT^ — -r± V x , Z0 .,y tZ0 \ = 0, (13) 

w • £ A,Z0 

which has the highest T c as solution for the critical temperature. This will be the physical solution. 

Results for an 11-layers system, obtained according to this scheme are shown in FIG. [| The critical tempera- 
ture T c is plotted as a function of the position of the nucleation point. The system parameters used apply for a 
V(100A)/Afc(100A) system, having a H C 2 t \\(0) of about 11 kGauss. The specific set of parameteres is given in the 
literatureflu 

For a small magnetic field, for 1 and 2 kGauss, the superconductivity occurs in the middle of the multilayer, as the 
highest T c is found for the nucleation point right in the center of the system. For such a small field the size of the 

system is of the same order of magnitude as £ = yj 2§jf i which quantity gives the scaling of the order parameter. For 1 

kGauss £ = 574 A, which is relatively large compared to the layer thickness. Microscopically, at low parallel magnetic 
field the pair amplitude is spread over the entire system and looks the same for any choice of the nucleation point 
zq. One could say as well, that for such small fields the entire system is just a surface. At higher parallel magnetic 
fields, the pairs are confined in about one layer around the position of zq. This makes that for higher fields a real 
difference can be expected for two different zq choices. For H > 3 kGauss, the nucleation point shifts from the center 
of the system towards the surface of the sample. For high fields one clearly sees, that superconductivity is found at 
the surface of the system only. 

The shift of the nucleation point towards the surface with increasing magnetic field illustrates the scale of the 
surface nucleation, which is given by the magnetic coherence length. For a bulk superconductor, if the magnetie-iield 
is applied parallel to the interface with the vacuum, the pair function is localized around zo nearby the surfaced 
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FIG. 6. The critical temperature versus the nucleation point, measured from the center of the layer. The magnetic field 
varies from 1 to 10 kGauss. Due to the symmetry of the multilayer, we only show one half of the system. 




FIG. 7. The points in which the nucleation occurs for different magnetic fields (with dots) fitted in terms of the magnetic 
coherence length £ (continuous line). 



F Se sp[-i(l_*)2], (14) 

in which z is now measured from the surface, situated at z = 0. The expression ( |l4| ) is a gaussian trial function which 
shows that the pair function is localized around the nucleation point zq . A detailed calculation which is based on the 
exact Weber-function solution shows that zq ~ £ and the optimal value is zq — 0.59 £. 

For our multilayer system, the dots in FIG. [t] show the nucleation point as a function of the magnetic field, which 
can be converted in a dependence on the magnetic coherence length. Measuring zq from the surface we find zq = C£, 
where C = 0.91, which is represented in FIG. ^ by the continuous line. It is not too surprising that this constant for 
the layered system is larger than the value of 0.59 for a bulk superconductor. The surface phenomena in the present 
systems will be shifted inward due to the covering N layer. 

The results presented up to now in this section apply to a multilayer system whose outer N layers have a thickness 
of lOOA, which is equal to the thickness of the S and N internal layers. As we were looking at the influence of the 
boundary on the superconductivity, we studied also systems in which the thickness of the outer N layers is larger. 
Keeping all other parameters constant, we only extend the thickness of the exterior layers to d™ ter = 350 A, such that 
it becomes larger than the magnetic coherence length £. The corresponding plot of T c versus zq is given in FIG. || for 
seven magnitudes of the magnetic field. Clearly, up to 5 kGauss the nucleation occurs in the middle of the system. 
No surface superconductivity does appear anymore since the normal outer layer is thick enough to "screen" it. For 
larger fields some surface superconductivity shows up again. 

In conclusion, in spite of the complexity of the layered systems, relations found between the nucleation and the 
scaling length £ of the superconductivity compare well with established results for homogeneous superconductors. 
This can be regarded as a support for the validity of the theory of Takahashi and Tachiki, which is devised to describe 
arbitrary inhomogeneous systems. 
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FIG. 8. For a thick exterior N layer, we have bulk superconductivity, rather than surface superconductivity. 

V. THE INFLUENCE OF THE PARAMETERS ON THE PROXIMITY EFFECT 

Up to now we have shown properties of finite multilayer systems by varying the applied magnetic field only. In 
order to prepare a solution to the as yet unsolved problem related to the position of the dimensional crossover in 
the H c 2 ,\\{T) curves, we need to study the way in which nucleation occurs by varying system parameters as well. 
In addition to the system parameters mentioned already, such as the diffusion constant far S and N layers and the 
density of states (DOS) of the S layer, we will consider the resistivity of the interfaces. Ejo Although we start from 
parameters which are specific to V/Ag, by varying some of them the set we are using might not characterise this 
system anymore. On the other hand, we are studying general properties of S/N systems, which means that we have 
to look not only to a particular system, but to a large class of S/N multilayers. 
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FIG. 9. The nucleation curves for a 11-layers S(120A)/N(120A) system, with increasing the S layer DOS. H=6 kGauss and 
Dn = 5 cm 2 /s. 

First we vary the S layer DOS in a S(120A)/N(120A) 11-layers system. This is a parameter which measures the 
density of electrons at the Fermi level. The higher it is, the more Cooper pairs can be formed, so the material is a 
better superconductor. While varying the S layer DOS, the other parameters keep the same values used before. In 
FIG. |^ we show the nucleation curves for an increasing DOS of the superconducting layer. The quantity dosS gives 
the ratio of the S layer DOS with respect to the N layer DOS, which covers a large range of DOS values, because in 
the equations only this ratio contributes. Since the system is symmetric with respect to the center of the multilayer, 
we only show half of it. First, we notice an increase of the multilayer critical temperature with increase of the S layer 
DOS, which is a consistent result. Secondly, with the increase of the S layer DOS, we notice a smooth transition 
from a regime which is dominated by the surface superconductivity, at lower dosS values, to a regime in which a 
modulation develops following the structure of the multilayer. Besides, surface superconductivity is still present by 
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a slight bump close to the surface. This bump cannot be present in infinite multilayers, having no surface. This is 
illustrated in FIG. |o[ in which the dosS = 25 result for the finite system, represented by the solid line, is compared 
to the periodic result for an infinite multilayer. The modulations are such, that the maxima are located at the centers 
of the superconducting layers. Apparently, for high dosS values the proximity effect is not strong enough anymore to 
smear out the high pair amplitude inside the S layers. This can be understood from the fact that at high dosS, the 
S/N interfaces behave like S/ Vacuum interfaces, so that the coupling between the S-layers is very small. 




FIG. 10. Nucleation curves for systems which have the same set of parameters, 
periodical multilayer, while the continuous one to a finite 11 layer system. 



The dashed line corresponds to the infinite 



Now we want to see how the behaviour displayed in FIG. [)] is modified by varying the diffusion constant Dm of the 
N layer and by looking at different magnitudes of the magnetic field. We concentrate on two representative curves, 
namely for dosS = 5 and 25. Results are shown in FIG. O. 

In FIG. 0(a) we show the evolution of the nucleation curve for dosS = 5 with increasing diffusion coefficient D^, at 
a lower magnetic field (H=2 kGauss) . Even the bump in the dosS = 5 curve of FIG. ^ has disappeared, because for such 
a weak field £ is large, and there is no surface superconductivity. For larger Dn values some surface superconductivity 
shows up again. This is because, in addition to the magnetic coherence length, from now on to be denoted by £ mas , 

the BCS coherence length comes into play. For dirty alloys this quantity is defined as £bcs — \/ JnkT ^ anc ^ ^ or ^ N 

= 20 cm 2 /s and T = 2 K it is equal to 348 A. This implies a good coupling between the superconducting layers, by 
which the entire multilayer develops properties pointing in the direction of a homogeneous superconductor. Because 
£mag = 406 A, for such a system surface superconductivity could develop at 0.59 x 406 = 239 A from the surface. 
The weak maximum in the D n = 20 curve in FIG. |ll](a) lies somewhat more inward. 

In FIG. |ll|(b) the magnetic field has a higher value, H=5 kGauss, which lies close to the value used in FIG. [|. Again 
dosS = 5, and the curve for Dn = 6 is similar to the dosS = 5 curve of FIG. ||. The effect of increasing Dn is much 
more pronounced than for the weaker field, and the maxima of the curves lie at a distance of about 0.91£ maff = 234 A 
from the surface. For large Dn values surface superconductivity is dominantly present. 

FIG. 0(c) shows the modifications to the dosS = 25 curve of FIG. [|. We restricted the calculation to a 7-layer 
system, by which one oscillation has dropped out. The magnetic field is H — 5 kGauss, slightly weaker than the field 
used in FIG. |[ The modulation disappears with increasing diffusion coefficient, because a higher D^ facilitates the 
coupling between the S-layers. 

All three panels of FIG. 11 show a decrease of T c with an increasing diffusion coefficient. This can be understood 
globally in terms of superconducting bulk properties. For a bulk superconductor the only eigenvalue that contributes 
in the determinant (^) is the ground state eigenvalue, which is proportioaaLtcuthe diffusion constant and the magnetic 
field, and inversely proportional to the magnetic flux quantum 4>q = ^,EJOEj 



2e ' 
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Following de Gcnnes we write cq in Eq. (Q) explicitly in terms of the so-called universal function Un(x) asEl 

'T C {HY 



e G = A:T c (0)Un 



Tc(0) 



(15) 



(16) 



while Un(x) is monotonically increasing with decreasing x. This implies a decrease of T C (H) with increasing ec, and 
therefore, according to Eq. (p"5|), with increasing D. Also, for a larger magnetic field the effect will be stonger. For 
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the layered systems an average value of Ds and Dm has to be used, but this average increases with increasing Dm- 
A decrease of T c with increasing Dm is seen indeed, and also the effect is stronger for a larger field. 

The scaling change due to the variation of £, ma g through its dependence on H is shown in FIG. fijl. The H = 6 
kGauss curve reproduces the dosS = 25 curve of FIG. ||. The nucleation curves show a smooth transition from curves 
which exhibit only surface superconductivity, to curves in which a modulation over the S and N layers is present. 
This transition happens when H« 5 kGauss and £, mag ~ 257 A. This value is comparable with the periodicity of the 
multilayer, d-^ + d$ = 240 A. In other words, when £ mag becomes as small as the periodicity of the multilayer, the 
nucleation exhibits the modulation corresponding to the geometrical structure of the multilayer. The maxima are 
located in the S layers and a preference is shown for the superconductivity to nucleate there. Comparing with Fig. 
^, in which this modulation appears when the coupling between the S-layers becomes weak at high dosS, here this 
feature of the nucleation curves can again be understood by the fact that the S-layer coupling decreases with the 
decrease of the coherence length at high magnetic fields. 

Finally we vary the boundary resistivity pb at the interfaces, expressed in units pQ.cm 2 . Since we have in view 
interfaces, we restrict our .calculations to just a thin superconductor film, covered by a normal metal. We vary ps 
around experimental dataJIS which are specific to Nb/Pd systems. We did calculations for a thin film of Nb(600A) 
covered by Pd(100A) at an intermediate magnetic field, which for Nb/Pd is about 12.5 kGauss. The nucleation curves 
are shown in FIG. |l3|. We notice that with increasing boundary resistivity ps, which means that the S/N interfaces 
become less transparent to the electrons, the critical temperature increases. This can be understood as follows. The 
interface resistivity makes the spreading of the electrons (Cooper pairs) over the different layers more difficult, which 
effectively leads to an enhancement of the pair amplitude in the superconducting layers. This weaker proximity effect 
results in a higher T c of the system. This effect is similar to what was discussed in Sec. [II regarding H&u, which 
systematically turned out to be larger than H C 2.±. 



VI. CONCLUSIONS 

The theory of Takahashi and Tachiki, applied so far either to a bilaye or to infinite multilayers,!^ has been 
applied to describe superconducting properties of finite multilayers. We have calculated the spatial dependent pair 
amplitude for two magnetic field directions, parallel and perpendicular to the layers. These properties, including the 
upper critical field curves, have been studied as a function of the number of layers. 

The nucleation problem, which appears when one studies systems to which a parallel magnetic field is applied, 
is illustrated on model systems such as V/Ag multilayers and Nb/Pd thin films for a large range of parameters. 
Moreover, we compare the results obtained for finite systems to what is found for infinite periodic multilayers. Since 
we take into account the boundaries, typical for a finite sample, we get surface superconductivity. 

Varying the parameters we see how they influence the proximity effect and the scaling of the sample and layer 
thicknesses to the magnetic and BCS coherence lengths. The superconductivity occurs differently for different choices 
of the parameters. If we define the dimensional crossover as being the temperature where the modulation over the 
geometrical structure of the multilayer appears in the nucleation curves, then we must expect a shift of the dimensional 
crossover in the H c2t \\ (T c ) curves under a change of the parameters. 

The results we get are meant to prepare an improvement upon the previous approaches in explaining the experi- 
mental results.0 In addition, the results obtained contribute to new confidence in the Takahashi- Tachiki theory. This 
comes from the observation, that our results about surface superconductivity in finite multilayers, being typically 
anisotropic systems, can be understood as a generalization of what is known about the phenomenon in isotropic 
homogeneous superconductors. 
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FIG. 11. The temperature as a function of the nucleation point zo, with increasing diffusion coefficient of the normal metal 
D N . (a)at H=2 kGauss, ^gg| = 5; (b)H=5 kGauss, §^§§§ = 5; (c)H=5 kGauss, = 25. 
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FIG. 12. The temperature curves for different nucleations, with increasing the external applied magnetic field H. 
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FIG. 13. The nucleation curves for different nucleations, with increasing boundary resistivity /)B(/ificm 2 ). The system 
thin film of Nb covered with Pd. The magnetic field applied is H=12.5 kGauss. 
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